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We summarize recent developments of the semiclassical description of shell effects in 
finite fermion systems with explicit inclusion of spin degrees of freedom, in particluar 
in the presence of spin-orbit interactions. We present a new approach that makes use 
of spin coherent states and a correspondingly enlarged classical phase space. Taking 
suitable limits, we can recover some of the earlier approaches. Applications to some 
model systems are presented. 

1. Introduction 

The periodic orbit theory (POT), initiated over 30 years ago by Gutzwiller, 1 is 
a semiclassical approach in which the level density of a quantum system is ap- 
proximated in terms of the periodic orbits of the corresponding classical system 
through the so-called 'trace formula'. It has provided a great stimulus to the re- 
search area of quantum chaos, ' 3 ' but is also applicable to integrable and nearly 
integrable systems. 5,6 ' 7,8 Although originally developed to describe the motion of 
a particle bound in a given external potential or an ideally reflecting boundary (a 
so-called 'quantum billiard'), the POT can also be applied to describe the quan- 
tum oscillations (shell-corrections) in many-fermion systems within the mean-field 
approximation. 6 ' 9 For a general introduction to the POT and its applications in 
nuclear, mesoscopic and nanostructure physics, we refer to a recent text book. 9 

An early application to nuclear physics consisted in a successful semiclassi- 
cal explanation of the systematics of ground-state deformations. 10 (A correspond- 
ing investigation was recently carried out also for metal clusters. 11 ) The onset 
of the mass asymmetry in the fission of actinide nuclei could also be explained 
semiclassically. 12 ' 13 Fig. 1 shows the shell-correction energy around the outer fission 
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barrier of 240 Pu, plotted versus the elongation parameter c and the mass asymme- 
try parameter a. The 3d plot to the left and the contour plots to the extreme right 
were obtained within the POT in terms of the (few) shortest periodic orbits, mod- 
eling the nucleus by axially symmetric cavities (see shapes to the extreme left), 12 
whereas the contour plots next to the right are the old quantum-mechanical results 
using realistic deformed nuclear shell- model potentials. 14 




Fig. 1. Shell-correction energy of 240 Pu near fission isomer and outer barrier vs. elongation c and 
mass asymmetry a. Left: 3d plot of semiclassical POT calculation, 1 along with the shapes at 
isomer minimum (A), asymmetric saddle (B), and on the way towards scission (C) (perpendicular 
lines indicate planes containing the shortest periodic orbits). Right: contour plot for neck parameter 
h = (upper panels) and h = —0.075 (lower panels). The left panels show the quantum-mechanical 
calculations with realistic Woods-Saxon potentials including pairing correlations and spin-orbit 
interaction; 14 the right panels show the semiclassical POT calculation for simple cavities. 12 



We can see that, in spite of the simplicity of the cavity model (with only one kind 
of nucleons, neglecting pairing and spin-orbit interactions), the semiclassical calcu- 
lation reproduces almost quantitatively the correct topology of the deformation 
energy surface and, in particular, predicts correctly the mass-asymmetric adiabatic 
fission path. The latter, indicated by the white dotted line, is simply given by the 
principle of stationary action of the shortest periodic orbits. The wavefunctions of 
the single-particle states, which are quantum-mechanically responsible for the low- 
ering of the barrier due to the asymmetry, were found to have their maxima exactly 
in the planes containing the shortest periodic orbits that dominate the asymmetry 
effect semiclassically. 13 

That the classical motion of the nucleons at these deformations is almost chaotic 
can be seen in Fig. 2. Here we show Poincare surfaces of section, obtained for 
angular momentum L z = 0, taken at the symmetric outer barrier (left) and at the 
asymmetric saddle (right); the corresponding shapes with the planes containing 
the shortest periodic orbits are shown on top of each plot. The regular motion is 
confined to small islands, containing the relevant periodic orbits, surrounded by a 
chaotic sea. Note that the energetically favored asymmetric shape has much smaller 
regular islands than the symmetric shape. Nevertheless, the shell effect coming from 
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these small regular regions of the classical dynamics is sufficiently strong to cause 
the collective asymmetry effect of the system. 




Fig. 2. Poincarc surfaces of section p (parallel momentum at reflection point) vs. 4> (polar angle 
of reflection point) for angular momentum L z = of the cavities corresponding to the symmetric 
outer barrier (left, c = 1.53, a = 0) and to the asymmetric saddle (right, c = 1.53, a = 0.13). 
The corresponding shapes and the planes containing the shortest orbits (vertical lines) are shown 
above the corresponding Poincarc plots. 13 

This example illustrates the strength of the semiclassical theory in explaining 
qualitatively (and, at times, even semi-quantitatively) important quantum effects, 
both in one-body and in many-body systems, in terms of classical dynamics. 

One aspect has, however, been neglected so far in the applications of the POT 
to nuclei: the spin of the nucleons. We know well that the spin-orbit interaction dra- 
matically modifies the shell effects that is, after all, why it had to be introduced to 
make the nuclear shell model work in the first place. 15 In the semiclassical calcula- 
tions of the Refs. 10 ' 12 ' 13 the neglect of the spin-orbit interaction was compensated 
by a simple readjustment of the Fermi energy, allowing one to locally reproduce 
the correct shell situations (or magic numbers). But this was only a temporary 
remedy, and a more rigorous semiclassical treatment of the spin-orbit interaction 
clearly remains highly desirable. The remainder of this paper is therefore devoted 
to some recent developments of the inclusion of spin degrees of freedom into the 
POT, which is by no means trivial since there is no classical analogon of the spin. 
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2. Semiclassical theories with spin 

We limit the discussion here to systems of fcrmions with spin s = 1/2 with a 
(mean-field) Hamiltonian linear in the spin operators 

H = H (r,p)+H so (1) 

with 

H (r,p) = H-+V(r), H so = hKvC(r,p). (2) 
2m 

The second term is a general spin-orbit interaction, where C(r,p) is an arbitrary 
vector function of coordinate r and momentum operators p, and k is a coupling 
strength independent of ft. a = (a x , o~ y ,o~ z ) are the Pauli matrices defining the spin 
operators s = ifter. In the non-relativistic reduction of the Dirac equation with an 
external electrostatic potential V(r), one obtains 

C(r,p) = [W(r) x p], k=1/4to 2 c 2 . (3) 

With the Coulomb potential this yields the familiar Thomas term for H so . 

2.1. Earlier approaches 

We first summarize two earlier semiclassical approaches. 

(i) SCL: Littlejohn and Flynn 16 developed a semiclassical theory for multi- 
component systems, treating the spin matrices quantum mechanically while Wigncr 
transforming the matrix operator (1) to the classical phase space (r, p), keeping the 
leading terms in an ft expansion. Diagonalisation leads to a pair of Hamiltonians 

ff±(r,p) = ffo(r,p)±ftK|C(r,p)|, (s = 1/2) (4) 

where H (r, p) and C(r, p) are the Wigner transforms (to lowest order in ft) of the 
corresponding quantum operators. H± can be considered as two classical adiabatic 
Hamiltonians with opposite spin polarizations. Their two sets of periodic orbits 
must be superposed in the final trace formula. This approach is often referred to 
as the "strong coupling limit" (SCL), since it becomes valid in the formal limit 
k — ► oo and ft — > with hn kept finite. 16 ' 17 The SCL approach suffers, however, 
from the problem of mode conversion (MC): whenever C = at a given point in (or 
in a subspace of) phase space, the two Hamiltonians H± become degenerate and 
singularities arise both in the classical equations of motion and in the calculation 
of the stabilities of the periodic orbits. (A similar situation occurs in the chemistry 
of molecular reactions when two or more adiabatic surfaces intersect.) The MC 
poses a difficult problem in semiclassical physics and chemistry that has not been 
satisfactorily solved so far. It was further discussed within the SCL in Rcfs. 18 ' 19 

(ii) WCL: Bolte and Keppeler 17 have derived a semiclassical theory from the 
Dirac equation. In the "weak coupling limit" (WCL) they arrive at a trace for- 
mula, in which the periodic orbits are given by the dynamics of the unperturbed 
Hamiltonian H and the effect of spin precession around the local 'magnetic field' 
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k C(r, p) appears through a simple modulation factor. This approach neglects terms 
of higher than first order in Kk and therefore is valid in the limit of weak spin-orbit 
couplings. Furthermore, it yields only a trivial spin degeracy factor of two in the 
trace formula whenever all periodic orbits of H$ are self retracing (i.e., librating 
between two turning points) such as in the two systems illustrated below. 

An instructive example is that of a two-dimensional electron gas with a Rashba 
type 20 spin-orbit interaction C = (—Tr y , tt x , 0) in an external homogeneous magnetic 
field where tt = p — eA/c. For this system the exact quantum spectrum is explicitly 
known, 20 and analytical trace formulae have been given for both the exact quantum- 
mechanical level density and the semiclassical WC and SC limits, 19 from which the 
limitiations of these two approches become evident. 

In a successful application of the SCL, a model relevant for nuclear physics was 
investigated in Ref. 19 It consists of a three-dimensional harmonic oscillator with 
Thomas-type spin-orbit interaction (3) 

V{v)= ]T ^ 2 r 2 , i = k<r.[VV(r)xp], (5) 

i—x : y : z 

which defines a Nilsson type 21 Hamiltonian appropriate for light nuclei (where the 
£ 2 term can be neglected). We express the oscillator frequencies in terms of two 
deformation parameters a, (3: 

u! x =oj , u) y = (l + a) ui , lu z = (1 + a) fj uj , (6) 

and use huj a s energy unit. For the general case of incommensurable frequencies 
(i.e., three-axial deformations), all periodic orbits of the unperturbed system Hq 
are librations along the coordinate axes and the WCL yields only the trivial spin 
factor of two. The shortest periodic orbits of the Hamiltonians H± in the SCL lie 
in the three coordinate planes and can be obtained analytically. 19 For these orbits 
no mode conversion takes place, and the SCL can be used. 

In Fig. 3 we show the shell correction to the level density of this system, with 
a deformation a = 0.1212, (3 — 2 and a spin-orbit strength k = 0.1 lo^ 1 . Both 
the quantum mechanical and the semiclassical 5g{E) have been coarse grained by 
Gaussian convolution over an energy range 7 = 0.5?iu;o, in order to suppress the 
contribution of the longer orbits and hereby to emphasize the gross-shell structure. 
The quantum-mechanical result is shown by the solid lines (QM) and includes the 
spin-orbit interaction in both curves a) and b). The semiclassical SCL result (SC) is 
shown by the dashed lines; in a) without spin-orbit interaction, which demonstrates 
that the latter dramatically changes the level density, and in b) with spin-orbit 
interaction. Only the six primitive planar orbits have been used. We see that this 
already leads to an excellent agreement with quantum mechanics, except at very 
low energies where semiclassics usually cannot be expected to work. As dicussed in 
Ref. 19 , bifurcations of the planar orbits occur for other deformations and values of 
k. These can, in principle, be handled by suitable uniform approximations (see Ref. 9 
for an overview), but they complicate the semiclassical calculations numerically. 
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Fig. 3. Coarse-grained level density Sg(E) of the three-dimensional harmonic oscillator (5) (see text 
for parameters, energy units: Itwo). Top and center: solid lines give quantum-mechanical results 
(QM) with spin-orbit interaction. Dashed lines (SC) give scmiclassical results (SC) obtained a) 
without and b) with spin-orbit interaction. Bottom: c) SC as in b) over a larger energy region. 19 

2.2. POT in an extended phase space 

A new semiclassical approach has recently been presented 22 ' 23 in which the spin 
degrees of freedom were introduced through spin coherent states 24 \z;s) defined by 

\z;s) = (1 + \z\ 2 )' s exp(zs + /ft) \s, -s) , s_|s, -s)=0, s± = s x ± is y , (7) 

where z = u — iv is a complex number. This allows one to define classical spin com- 
ponents n = (n X7 n y ,n z ) — (z; s\s\z; s) jhs and to enlarge the classical phase space 
by only one pair of canonical variables (w, v) , independently of the value of the spin 
s. Starting from the path integral in the SU(2) spin coherent state representation 25 
and making the usual stationary-phase approximation in its evaluation, the semi- 
classical dynamics of the system in the extended phase space (EPS) (r, p, v, u) is 
then determined 22 by the Hamiltonian 

H(r,p,v,u) = H (r,p) + h~K2sn(v,u) ■ C(r,p) . (8) 

Solving the equations of motion following from (8), one can determine the peri- 
odic orbits in the EPS and their properties, yielding the required input into the 
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Gutzwiller trace formula. 1 Special attention is required for the Maslov indices 26 
and other phases arising in connection with the spin degrees of freedom. 23 ' 27 

The EPS approach is free of the problem of mode conversion and, in principle, 
applicable for both weak and strong spin-orbit interactions. Note that the Hamil- 
tonian (8) explicitly couples the orbital degrees of freedom (r, p) with the spin 
degrees of freedom (u, v). This is illustrated for the following model Hamiltonian of 
a two-dimensional semiconductor quantum dot with Rashba interaction: 

H=(P 2 X + P 2 y )/2m* + m*(cj 2 x x 2 + iJ 2 y y 2 )/2 + Hk {a yPx - a x p y ) , (9) 

where m* is the effective mass of the conduction electrons, and the deformation 
parameters were chosen as lo x — 1.56 u>o and u> y — 1.23 u>q. (In the figures, units are 
chosen such that h — m* = ui n = 1 and E and k become dimensionless.) Like in 
the three-dimensional system discussed above, the periodic orbits of Hq are pure 
librations, so that the WCL yields only trivial results. But here also the SCL cannot 
be used, due to the MC problem, so that a new treatment is required. 

For < k < 0.7, the following set of 12 shortest periodic orbits in the EPS were 
found: 22,23 (i) Two pairs of orbits A^ and librating along the x and y axes with 
fully polarized spin n y = ±1 and n x = ±1, respectively, (ii) Two pairs of orbits 
and D x2 oscillating around A^ with n y ~ 0, and two pairs of orbits and 
Dy 2 oscillating around A y with n x ~ 0. For stronger couplings with k > 0.7, new 
orbits bifurcate from the A orbits. 

In Fig. 4 we show the (x, y) shapes of the orbits A+, D^, and D+ 2 (left panels), 
and the time dependence of their spin components n x , n y , and n z over one period 
(right panels), all evaluated for n — 0.67 and E = 60. We see that along the orbits 
D xl and D^ 2 , the spin rotates mainly near the (n X7 n z ) plane (i.e., n y ~ 0), but in 
a non-uniform way. This complicated spin motion, together with the wiggly orbital 
shapes of the D orbits, reveals the rather sophisticated dynamics which is obtained 
from the equations of motion in the EPS through the explicit coupling of spin and 
orbital degrees of freedom. 

Fig. 5 shows the shell correction of the level density, obtained quantum-mech- 
anically (solid lines) and semiclassically (dashed lines), both coarse-grained over an 
energy range 7 = 0.6. The upper panel contains the semiclassical EPS result, while 
the lower panel exhibits the WCL result which is identical to that obtained by ig- 
noring the spin-orbit interaction and multiplying the level density by a spin factor 
of two. (Note that the semiclassical trace formula for the unperturbed harmonic 
oscillator is analytically known and quantum-mechanically exact. 9 ) We observe a 
reasonably good agreement of the EPS result with quantum mechanics. The semi- 
classical amplitudes are too large, which is attributed to close-lying bifurcations of 
the periodic orbits; 23 suitable uniform approximations are expected to remedy this 
defect. The phases of the quantum oscillations are, however, very well reproduced, 
which is not achieved at all in the WCL. 
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Fig. 4. Periodic orbits in the two-dimensional quantum dot (9) (see text for parameters). Left 
panels: orbits in the (x,y) plane. Right panels: Spin components n x , n y , and n z versus time. From 
top to bottom: orbit A£ along x axis with polarized spin in y direction, and orbits D+ 2 and D^" 2 
oscillating around the x axis with spin rotating near the (n x ,n z ) plane. 22 ' 23 




Fig. 5. Coarse-grained level density of quantum dot (9) for the same parameters as in Fig. 4. Solid 
lines: quantum-mechanical result. Dashed lines: semiclassical results; in upper panel: EPS result, 
in lower panel: WCL result. 23 
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3. Summary and Outlook 

After a short review of the two earlier semiclassical approaches including spin- 
orbit interactions, corresponding to the weak coupling limit (WCL) and the strong 
coupling limit (SCL), we have presented a new approach that makes use of spin 
coherent states and leads to (semi)classical dynamics in an extended phase space 
(EPS). Both the WCL and the SCL (with the 'no- name' phase still lacking) could 
be recovered from the EPS approach taking suitable limits. 

We have only discussed here Hamiltonians linear in spin and systems with spin 
s = 1/2. The EPS approach was formulated for arbitrary spin-dependent Hamil- 
tonians and arbitrary values of s. As is well known, semiclassical approximations 
work best in the limit of large quantum numbers. Whether they can be used for 
small quantum numbers is a matter of numerical experience and fortune (like for 
harmonic oscillators without spin) . As to the spin, there are firm grounds to expect 
good semiclassical results also for s = 1/2 for Hamiltonians linear in spin. 23 

The example shown in Figs. 4 and 5 suggests that the EPS approach has a good 
potential for a general semiclassical theory with spin. Practically, this approach 
suffers from a large number of bifurcations occurring for the periodic orbits under 
variations of both k, energy, and deformation parameters. It would therefore be 
desirable to use it as a formal starting point for further approximations. 

Indeed, in the limit of a weak spin-orbit coupling such that H so <C H 0} the 
WCL trace formula of Ref. 17 could be rigorously derived from the EPS approach. 28 
Hereby the equations of motion of the orbital and spin degrees of freedom were 
decoupled, and the effects of H so were included in the phases of the trace formula 
in first-order perturbation theory. 29 Pushing this treatment to second order, along 
the lines proposed in Ref. 30 for spin-independent systems, might allow one to extend 
the WCL trace formula, so that it becomes valid for larger values of k while still 
benefitting from a simpler determination of the periodic orbits. In particular, for the 
self-retracing orbits where the first-order result is trivial, the second-order treatment 
is expected 30 only to lead to a phase correction while the amplitudes in the trace 
formula still are determined by the unperturbed orbits of Hq. This expectation is, 
indeed, strongly supported by the results shown in the lower part of Fig. 5. 

On the other hand, a careful study 23 of the situation H so ~ Hq reveals that 
the essential ingredients of the SCL approach of Refs. 16 ' 18 - without, however, the 
so-called 'no-name phase' - can also be retrieved from the EPS approach. Adding 
suitable corrections to this limit might help to overcome the mode conversion prob- 
lem without going through the cumbersome task of finding all relevant periodic 
orbits in the extended phase space, including all possible bifurcations and their 
treatment by uniform approximations. 

In a recent application of the EPS approach to mesoscopic transport theory, the 
effects of spin-orbit interaction on weak anti-localization have been investigated, 31 
and its application to semiclassical studies of nuclear shell structure including re- 
alistic shell-model potentials and spin-orbit interactions is in progress. 32 
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